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There is disagreement in the literature as to how the center of flexure of an isotropic 
elastic beam should be denned. In either of the two most widely accepted definitions, upper 
and lower bounds for the coordinates of the center of flexure are obtained. 

1. Introduction 

We consider a cantilever beam of uniform cross section and choose the 2-axis to coincide 
with the line of cross-section centroids. The plane end z=0 is assumed fixed and the other 
end z=l is loaded by some distribution of forces which is statically equivalent to a single force 
with components (W x , W y , 0) acting at a load point (x L , y L , /). In the treatment of the flexure 
problem [9] 1 (a similar treatment appears in [3]) we seek a load point (x r , y f , I) which has the 
property that any load (W x , W y , 0) applied through this point produces no local twist at the 
centroid of the section, i.e., the mean value of the local twist over the cross section vanishes 
(see [9] p. 200). This point (x f , y f , /) is referred to as the center of flexure [9]. The determina- 
tion of (x f , y f , I) permits us to break up the general flexure problem into two separate problems — 
one of pure torsion and one of pure flexure. This point is sometimes referred to as the center 
of shear (see, for instance, [3]), and it is well known that there is difference of opinion in the 
literature as to the way in which the point should be defined (see, for instance, [1, 2, 4, 10, 11]. 
A discussion of these differences has been given in a recent paper by Pearson [7], who demon- 
strates that in general the two expressions for the center of shear yield results which are very 
near to one another. 

It is not our point here to enter into the discussion of the relative merits of the two defini- 
tions. We merely point out that the expression of Trefftz [10] is simpler and independent of 
Poisson's ratio. The center of flexure of [9], on the other hand, permits an easy decomposition 
of the general flexure problem. 

In order to avoid confusion we shall hereafter refer to the point defined in [9] as the center 
of flexure and that given by Trefftz [10] as the center of shear. 

In this paper then we derive upper and lower bounds for the coordinates of the center of 
flexure (x f , y f , I) and the center of shear (x s , y s , I). In the first case we use the formulation of the 
flexure problem derived in [9] and thus concern ourselves with the problem of pure flexure. 

It is of course well known (see, for instance, [7]) that if we have the solution to the pure 
torsion problem for the beam in question then we can determine both the center of shear and 
the center of flexure without knowledge of the solution to the flexure problem. As would be 
expected it is possible in each case to obtain upper and lower bounds for the coordinates of the 
point by approximating the solution to the torsion problem alone. 

2. Bounds for the Center of Flexure 

We denote by & the cross section of the beam and let C be its boundary. (The region 
may be simply or multiply connected.) The origin of a rectangular coordinate system is 
taken at the centroid of @. For simplicity the x and y axes are chosen to be principal axes of 
inertia, i.e., 



9 



xydA=0. ^ 



* Figures in brackets indicate the literature references at the end of this paper. 
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According to [9] the pure flexure problem is solved once the solutions to the following boundary 
value problems have been determined : 

*HS-0. &+&=» ^3 (2.2, 



d^ 2 d?/ 2 ' dx 2 d?/ : 



^g=[(l +ff )z 2 -<^K (2.3) 

^=[(l+^)2/ 2 -« 2 K on (7 

where a is Poisson's ratio and n x and t^ are the x and ?/ components of the unit normal on C. 
(The normal is assumed directed outward from @.) The coordinates of the center of flexure 
(x f , y f ) are defined by : 

*<-mh)z j> tr y ^'^^ (2 - 5) 

= f(W (2.6) 

I y = ffx 2 dA. (2.7) 



where 



:= 



^ 



Thus 



2(l+<r)i> / =JJ[|^ (x%)-^ (^-(l+^+^JcM. (2.8) 

Using the divergence theorem we have 

2(l+<r)I x T f =(D<p 2 (xn v -yn x )ds-J (2.9) 

c 
where 

J= \\[{\-\-a)xy i —<rx z \dA. (2.10) 

Since J is a geometric quantity it can be computed explicitly. 

We note now that the coefficient of <p 2 in the boundary integral in (2.9) is proportional 
to the boundary value of the normal derivative of the warping function <p in the pure torsion 
problem. In fact 

x—=yn x —xn y on 0. (2.11) 



Thus if we let 
we find 

where 



X=2(l+ff)I x x f +J (2.12) 

x =-f <P2*~ds=-D(<p,<p 2 ) (2.13) 



The determination of bounds for x r is thus equivalent to the determination of bounds for x- 
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Note that if the warping function <p is known the value of x f can be determined explicitly 
since 

X=-6<p[(l + a)y 2 -ax 2 }n y ds. (2. ir>) 

c 

If <p is not known we introduce two arbitrary harmonic functions u and u 2 and define 

&!= — D(<p,u 2 ) = — CD u 2 (yn z — xn^ds, (2 . 1 G) 

c 

a 2 =—D(u,<p2) = — (bu[(l + a)y 2 —a2?]njls, (2.17) 

c 

a s =-D(u,u 2 ). (2.18) 

With u and u 2 prescribed we can determine a u a 2 , and a 3 explicitly. We now form the expression 

X—a l —a 2 +a 3 =—D((p—u,<p 2 —u 2 ) 

= -§(*-«) «f [(1+ ffV- <rZ 2 K~ ~ 2 \ A. (2.19) 

An application of Schwarz's inequality yields 

c c ' 

It remains now to approximate the first boundary integral on the right of (2.20) in terms of 
known quantities. 

By adding an appropriate constant (if necessary) it is possible to choose <f-u in such a 
a way that 



)(<p-u)ds=0. (2.21) 

c 

With ip- — % so chosen we shall estimate the first integral on the right of (2.20) in terms of the 
boundary integral of the square of its normal derivative. To this end we note that 



M <—>)'" f(M) 



2 

ds 



C 47T 2 



J__ ->min^— 5 =^ ( 2 - 22 ) 

(D(<p-u) 2 ds £*-o &)v 2 ds 
c c 

where L denotes the length of the perimeter of C, bjbs denotes the tangential derivative on 

C, and v is any function continuous on C and satisfying the condition (J) vds = 0. The quantity 

c 
on the right of (2.22) is just the eigenvalue of the vibrating string equation, the differential 

equation which arises as the Euler equation for the minimum of the quotient ( (D(()fl/d.s j ) 2 ds V 

c 
(yv 2 ds. It follows then that 
c 

(Jj(»,-„)%<g j> [A {<p - u )Jd S . (2.23) 

c c 

We seek now an inequality which relates the integral on the right-hand side of (2.23) to the 
boundary integral of the square of the normal derivative of <p—u. We introduce therefore, 
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an arbitrary continuous vector function with components j k and consider the generalized 
Rellich identity (see [6, 8]). Letting w=(p—u, we have 







3 



c c q q, (2>24 ) 



•S AlRLi "9 >y *SS* 



where A denotes the Laplace operator (A=(d 2 /c);r 2 ) + (d 2 /d?/ 2 )) and summation is to be 
carried out over the repeated indices &,Z=1,2. The f k are now assumed to be so chosen that 



p=Jkn k >0 on C. 



(2.25) 



A decomposition of the derivatives of u in the boundary integrals into normal and tangential 
components yields (see [8]) 

*c)f k ()w dw 



HW^-HW^f^^^M^^-^! 



ds ds dn 



$ 



bxi bx k bxi 



dA. 



(2.26) 



Since the functions j k are prescribed, we can find a constant y such that (see [6, p. 555]) 



Sf^ 



<ta* 



2 ffM 

J J dxidxjcdxt 

We make use also of the ordinary Green's identity 

D(w,w)=(j)w-^ds, 
c 
which yields by Schwarz's inequality 

D{w,w) < (jwhlsj^ji^jdsj- 
c c 

In view of (2.23), eq (2.29) yields 



(2.27) 



(2.28) 



(2.29) 



(2.30) 



We now insert (2.27) and (2.30) back into (2.26), make use of the Schwarz inequality and 
obtain 






yL 



2 m 



mm^mw ^ 



where p min is the minimum value of p on C. This inequality yields 

f /r 



f ^% 



tfpffi'lAfif^'H&fffi*! 



+< 



Q-W 






ds 
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From (2.32) it is thus always possible to obtain the inequality 



where 
and 






3=< 



ds 



yL 



max -TT/^min 



(2.33) 

(2.34) 

(2.35) 



If the boundary C is star-shaped with respect to some point (x ,y ) in @, the problem is 
simplified considerably. The fact that the boundary is star-shaped with respect to (x ,y ) is 
equivalent to the condition that 

<L=(x-x»)n x +(y-y,)n y >0 (2.36) 



at every point on C. Thus, in this case, (2.25) is satisfied if we let 



(2.37) 



The left-hand side of (2.27) then vanishes identically, and hence we take y = in (2.31), (2.32), 
and (2.35). We find in this case instead of (2.33) (see [5]) that 



C C 

In the general case then, an inequality of the form 



(2.38) 



(2.39) 



is obtained, where A is a completely determined positive function. Inserting (2.39) and 
(2.23) into (2.20), we find 

(x-a 1 -a 2 +a,y< gj ^ A g-^+^J \ls £ /^-[(l+^-rflnA^ (2.40) 
c c 

(It is apparent that the quantity A may be placed in the second integral rather than in the 
first if it is desirable.) We now use the Rayleigh-Ritz technique to make the two integrals 
on the right of (2.40) small. 

Note that it is possible to derive close upper and lower bounds for x f without approximating 
the flexure function at all. In fact, if we use instead of u a conjugate function v defined by 2 

du dv du dv , A 1 v 

x-=sr> ^= — ^-> (2.41) 

dx dy by dx 

we obtain (setting ^ 2 = 0) 

|x-a 2 | 2 < g ^[A («-|r 2 )] 'da §A{[(l-*)y*-<rx*}n v }hls. (2.42) 

c c 

In this case the A term has been put in the second integral since this expression is a purely 
geometrical quantity (independent of v) and can be computed explicitly. 



The function v is not to be confused with that used in (2.22). 
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In an entirely analogous way we find that 

-2(l + <r)I y y f +jj[(l + <r)x 2 y-cry*]dA=D(<p,<p 1 ). (2.43) 

By introducing the two arbitrary harmonic functions u and u x and proceeding as before r 
we obtain upper and lower bounds for y f . Note that if we choose %,%s0 and merely approxi- 
mate the warping function (or the conjugate function) on the right, then the same minimizing 
function u (or v) minimizes the error term in each case. 

3. Eounds for the Center of Shear 

We look now at the expression for the center of shear (x s ,y S) l) as defined by Trefftz [10]. 
The coordinates x s and y s are given by 



9 

Equation (3.2) may be rewritten as 



x s =-j j LydA (3.1) 

* 9 

j- \LxdA. (3.2) 



V 'SL 



jj- fj(Aa?) v dA. (3.3) 

* 9 

By Green's formula we have 

c 

=77-0) [3x 2 n x <p— x z (yn x — xn y )} ds. (3.4) 

y c 

The divergence theorem then yields 

Vs+2J I Wv dA =2j ty x 2 nx<pds. (3.5) 

We introduce the notation 

Q^Vs+^j jjtfydA, (3.6) 

and two harmonic functions h x and h 2 . Let 

b^&h^rizds, (3.7) 

c 

b 2 =(b<p^ ds=(bh2[yn x —xny\ds, (3.8) 

c c 

=fk^ds. (3.9) 



^3 

C 
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Combining (3.5) to (3.9) we obtain 

$-b l -b 2 +b. s =(b(<p-h l )(x 2 ri x -h 2 )ds. (3.10) 

c 

Then by Schwarz's inequality 

[*-J 1 -6 f -6,]»<^(^-A 1 )WB(f(a^,--^Y(fe. (3.11) 

c c 

Again it is desirable to choose (<p—hi) in such a way that 



f 



(<p—h 1 )ds=Q. (3.12) 



Using (2.23) and (2.39), we obtain finally 

[ »- S| - J , +W <g^g-^)'<^( A ,-^)'*. (3.13) 

c c 

We insert the known boundary value for <p on the right-hand side and use the Rayleigh-Ritz 
technique to obtain close upper and lower bounds for y s . Note that if the boundary value 
problem 

AA 2 =0 in <2i 

^=* 2 n x on C (3.14) 

can be solved, it is then possible to determine y s explicity without solving either the torsion or 
the flexure problem. 

In a similar way it is possible to obtain upper and lower bounds for the coordinate x s . 
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